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ABSTRACT 

We carry out an exploratory weak gravitational lensing analysis on a combined VLA 
and MERLIN radio data set: a deep (3.3/iJy beam -1 rms noise) 1.4 GHz image of the 
Hubble Deep Field North. We measure the shear estimator distribution at this radio 
sensitivity for the first time, finding a similar distribution to that of optical shear esti- 
mators for HST ACS data in this field. We examine the residual systematics in shear 
estimation for the radio data, and give cosmological constraints from radio-optical 
shear cross-correlation functions. We emphasize the utility of cross-correlating shear 
estimators from radio and optical data in order to reduce the impact of systematics. 
Unexpectedly we find no evidence of correlation between optical and radio intrinsic 
ellipticities of matched objects; this result improves the properties of optical-radio 
lensing cross-correlations. We explore the ellipticity distribution of the radio counter- 
parts to optical sources statistically, confirming the lack of correlation; as a result we 
suggest a connected statistical approach to radio shear measurements. 
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1 INTRODUCTION 

Gravitational lensing is a powerful probe of the large 
scale structure of the Universe, due to its ability 
to study the matter distribution without differentiat- 
ing betw een luminous, baryonic matte r and dark mat- 
ter (e.g. Bartelmann fc Schneiderl l200ll ; iRefregierl l2003bl; 



ISevmour et al.l [20041 ). leading to a comparable source den- 
sity to that found in the optical. An advantage of radio over 
optical is that radio interferometers have well known and 
deterministic beam patterns which may allow for more ac- 
curate modeling of the effective PSF. 

The only detection of c osmic shear at radio wavelengths 



iMunshi et al.ll200St ). The effect has the added benefit of be- 
ing sensitive to the geometry of the universe, makin g it use- 
ful for the study of the dark energy (|Hutererll2002h . 

Almost all cosmic shear analyses so far have been car- 
ried out using optical data. However, future radio surveys 
could be competitive for lensing when compared to their 
optical counterparts, for several reasons. Firstly, radio sen- 
sitivities are rapidly reaching a level at which the radio 
emission from ordinary galaxies will be routinely resolved 
(e.g. with e-MERLHSQ LOFAE0, and eventually SKA0; c.f. 
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so far was car ried out by Chang et alj (12004) , usi ng the 
FIRST survey (|Becker et al.lll995i ; IWhite et alill997l ). This 



has a detection threshold of 1 mjy, with a consequent ~ 20 
resolved sources per square degree useable for weak lensing; 
this is a much lower number density than found in deep op- 
tical shear surveys, with a concomitant lower signal-to-noise 
on the final cosmological constraints. However, the differ- 
ential radio source counts at 1. 4 GHz show an incre ase at 
flux densities below 1 mjy, e.g. ISevmou~ et all (|2004l ), and 
it is this increase in the number density at the micro- Jansky 
level that makes future radio weak lensing attractive. 

Given this, it is of interest to make an initial foray into 
deep radio weak lensing with some of the most sensitive 
radio data currently available. Using MERLIN and VLA 
deep observations of the Hubble Deep Field North (HDF-N) 
jMuxlow et al.ll2005l ). we seek to understand the properties 
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of the radio source ellipticity distribution down to ~ 50 fj,Jy, 
including systematic effects, and to make radio shear mea- 
surements in this field, providing cosmological constraints. 
Along the way, we will compare the radio shapes with corre- 
sponding shapes measured in the HDF-N optical image ob- 
served with the HST ACS (|Giavalisco et al.ll2004 ), in order 
to examine the power of cross-correlating shear estimators 
from the optical and the radio. This modest investigation 
will therefore pave the way for more ambitious radio weak 
shear surveys (and joint radio-optical surveys) in the near 
future. 

Our paper is organised as follows: in Sj2]we outline the 
theoretical background and notation for weak gravitational 
lensing which we will use throughout, and then introduce 
the data used in 331 In we describe the methods used 
to measure shear estimators for each source detected in our 
image; then in ^5]we present the results of this process, giv- 
ing shear estimator histograms for radio and optical data, 
and examining the impact of systematic effects. We present 
shear correlation functions for the radio and optical data, 
together with optical-radio shear cross-correlation functions, 
providing constraints on cosmological parameters. In ij6]we 
compare the optical and radio shear estimators, finding lit- 
tle correlation between the intrinsic ellipticities of objects in 
the two bands; we show how this adds to the usefulness of 
cross-correlating optical and radio shear estimators. In Sj7]we 
present an alternative statistical method of cross-correlating 
the two datasets, before drawing conclusions in iJH] 



2 WEAK GRAVITATIONAL LENSING 

In this section we briefly describe the relevant theory be- 
hind weak gravitational le nsing and the cosmic shear statis - 
tics used i n this w o rk (e.g. Bartelmann fc Schneider! (|200lf ) , 
iRefregieri (|2003bh . iMunshi et alj (|2008T ) for comprehensive 
reviews) . 

As light travels from a distant source to an observer 
it is deflected by the tidal gravitational field of the inter- 
vening matter. This deflection angle maps the position of 
the source j3 on the source plane to position 8 on the image 
plane. Since the deflection angle varies from place to place, 
extended sources like galaxies suffer a change in their ob- 
served shape; in the weak lensing limit this effect can be 
quantified by a 2 x 2 distortion matrix 
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where k is the convergence, which describes the isotropic 
change in the object size, and 71 and 72 are shear compo- 
nents describing anisotropic stretching of the object along 
the x-axis and at 45° to the a;-axis. di is the partial deriva- 
tive with respect to 8i, and <j> is the lensing potential, related 
to the Newtonian gravitational potential $ via the line of 
sight of sight integral 
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where \ i s the comoving distance. The shear and conver- 
gence can then be written as 



7i = - 9 2 )<t>, 72 = did 2 <t>- 



(3) 



For the statistical analysis of cosmic shear we require suit- 
able 2-point statistics. Here we only consider a 2-dimensional 
shear field (i.e. projected on th e sky), a l thoug h a full 3- 
dimensional field can be used, e.g. lHeavensl |2003l ). Following 
iMiralda-Escudd i|l99ll ) we define tangential and 45° rotated 
shear components, 7t and y r , in relation to the line connect- 
ing a pair of galaxies: 



7t 



71 cos 28 + 72 sin 28 y r = 72 cos 28 — 71 sin 28, (4) 



where 8 is the position angle between the a;-axis and the 
line that connects the two galaxies. There are three pair- 
wise shear correlation functions 
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where averaging is done over all pairs of galaxies separated 
by an angle 8 = \6\; in practice, we will average over pairs 
in a bin, i.e. pairs with separation 8 ± Ad. Due to parity 
symmetry £ tr = f rt should equal zero, and this can be used 
to test for residual systematics. We further define the corre- 
lation functions £±(0) which are related to the convergence 
power spectrum P K by 

£±(0) ee &($) ± £„.(<?) = J ^P*(£)JoA&) M ( 6 ) 

where J n is a Bessel function of order n. The convergence 
power spectrum is a projection o f the matter power spec- 
trum Pg along the line of sight e.g. iBartelmann fc Schneider] 
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where \ is the comoving distance along the line of sight and 
Xi is the maximum comoving distance of the survey, fx (x) is 
the comoving angular diameter distance, n(x) is the source 
distribution and £ is the wavenumber on the sky. 

Weak gravitational lensing arises from scalar pertur- 
bations of the space-time metric, and therefore the shear 
field is expected to possess no handedness. The statistics 
of this scalar field can thus be described by a single corre- 
lation function, meaning that there must exist a degener- 
acy between the two c orrelation functions in equation ©. 
ICrittenden et alj (|2002f ) have shown how to transform these 
two correlation functions into gradient (_B-mode) and curl 
(B-mode) components. The _B-mode contains the true lens- 
ing signal while the -B-mode is expected to be zero in the 
absence of noise. Clearly, a non-zero .B-mode indicates a 
non-gravitational contribution to the sh ear field and can 
theref ore b e used as a t est fo r systematics. ICrittenden et al.1 
|2002i ) and |Pen et all (|2002i 'l derive the decomposition for 
the shear correlation functions, 
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where the definition of £'(#) is also given in lSchneider et al] 
(|2002h 
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In practice data will only contain information out to some 
radius 61, and so the integrals will be truncated. 



3 DATA AND CATALOGUE CREATION 

In this section we describe the radio and optical observations 
we will analyse, and give details of how the source catalogues 
were prepared. 

3.1 Radio data 

The HDF-N has been subject to some of the deepest ra- 
dio observ ations made to dat e using both the VLA and 
MERLIN (jMuxlow et al.ll2005l L In this study we have used 
data combined from these two arrays to provide both high 
sensitivity and sub-arcsecond imaging across a large field- 
of-view. These comb ined data were originally presented by 
iMuxlow et al.l (|2005l ) . For completeness in the following sec- 
tion we include a summary of the basic characteristics of 
these observations, the data processing and creation of the 
science images used for this study. 

3.1.1 Radio Observations 

The VLA in A-array observed the HDF-N for a total of 
50 h at 1.4-GHz in November 1996. A detailed description 
of these observ ation s and their cal i bratio n can be found in 
iRichardsl (|200Ch and lMuxlow et all (120051 ). To provide wide- 
field imaging capabilities and minimize the effects of chro- 
matic aberration these data were observed using the VLA in 
'pseudo-continuum mode', correlating parallel hands of po- 
larisaton into 7x 3.125-MHz channels, centred on intermedi- 
ate frequencies of 1365 and 1435 MHz. These data were sub- 
sequently time averaged to 13 s and the full field to the first 
side-lobe of the primary beam was imaged at low angular 
resolution, employing significant spatial tapering, to identify 
bright, confusing sources which might contaminate the inner 
portion of the field. Each identified confusing source, above 
0.5 mjy, was then carefully imaged in full spectral-line mode 
to account for changes in the primary beam response over 
the passband and then subtracted from the uv data. Over 
the course of the 50 h observations several isolated time seg- 
ments were identified where it was not possible to adequately 
subtract these confusing sources. This incomplete subtrac- 
tion of confusing sources in individual time segments was di- 
agnosed as probably associated with time variable telescope 
point errors, or possibly low-level interference. All data as- 
sociated with these times was deleted, resulting in a total of 
42 h of high-quality data which was used for all subsequent 
imaging and combination with the MERLIN data. 

Between 1996 February and 1997 April the MERLIN ar- 
ray, including the 76-m Lovell telescope, was used to observe 
the s ame location for a total of 18.23 d (see IMuxlow et al.l 
120051 ). All observations were made at 1.40 GHz, with the par- 
allel hands of circular polarisation correlated into 32 chan- 
nels of width 0.5-MHz with 4 s time averaging. This config- 
uration allowed imaging over a 10- by- 10 arcmin field, com- 
parable to the half-power beamwidth of the Lovell telescope 



(12.4 arcmin). Due to the lack of short interferometer spac- 
ings in the MERLIN data and the primary beam of the 
Lovell telescope, confusing sources beyond the 10-arcmin 
field caused no significant problems. However to mitigate 
any potential problems the four brightest sources within the 
10-arcmin field were mapped and subtracted from the vis- 
ibility data prior to combining these MERLIN data with 
data from the VLA. 



3.1.2 Data Combination and Image Creation 

For both the MERLIN and the VLA data the observations 
and correlator configurations were optimised for each of the 
separate arrays; maximising the sensitivity and field-of-view 
available and minimising the effects of time and chromatic 
aberration. Consequently each of these two data-sets con- 
tains fundamentally different configurations including dif- 
ferent intermediate frequencies and channel configurations. 
Due to these data structure differences, a combination of 
computing and software limitations, and the scientific re- 
quirements to image a very large field-of-view the combina- 
tion of these two data-sets in the visibility plane was not 
feasible. 

Following and extending the method outlined in 
IMuxlow et all (|2005h these MERLIN and VLA data were 
imaged and combined in the sky plane completely covering 
a mosai c area of 8.5-by-8.5 arcmin centred on th e HDF-N 
(see also lRichards et al.| [2007; Bcs wick et al.|[2008h . Initially 
small individual image facets were created from both the 
VLA and MERLIN data at a series of positions completely 
covering the 8.5-by-8.5 arcmin field separately and with no 
deconvolution applied. In each case these individual, 'dirty' 
map and beam facets were created by naturally weighting 
the data and using the optimum griding spacing for each of 
the arrays (MERLIN cellsize 0705 and VLA 074). Follow- 
ing the re-griding of the coarse sampled VLA 'dirty' maps 
and beams, the individual 'dirty' map and beam facets of the 
MERLIN and VLA data were averaged to produce combined 
images with no deconvolution applied. The central quarter 
of these combined dirty map facets were then deconvolved 
with the combine d-array dirty b eams using a conventional 
CLEAN algorithm (jHogbom 1974]) as implemented within the 
AlPS routine APLCN. The resultant facets were restored using 
a circularly symmetric Gaussian beam of width 074. Each 
individual deconvolved facet was then mosaiced together, 
fully accounting for sky-curvature effects, to create a single 
clean image of the field. 

Extensive tests comparing this image based data combi- 
nation with a visibility plane combi nation of small e r mos aics 
from the same data are detailed in lMuxlow et al.l (|2005l . see 
their Figure 1). These tests show that these two approaches 
produce essentially identical results. 

3.1.3 Radio Catalogue 

For both radio and optical source e xtraction, we use 
the f reely available software SExtractor l|Bertin fc Arnout3 
Although SExtractor was devel oped for the anal- 
ysis of optical d at a, several authors (|Bondi et al.l 120031 ; 
iGarn et al]|2008al rbl; iHuvnh et all 120071 ) have shown that it 
is able to generate reliable noise maps and locate objects 



4 Patel et al. 



within radio images. It should be noted that SExtractor is 
highly sensitive to input parameters and there is no single 
output catalogue that will be suitable for all applications. 

We make two radio catalogues. Firstly we create a 
more conserva t ive 'g old' set based on the catalogue of 
iBiggs fc Ivisonl (|2006l ). That work used the 10' x 10' VLA 
pointing alone with an independent reduction, detecting ob- 
jects at 5cr in the VLA image alone. Their catalogue contains 
537 sources, 83 of which are contained within the co-added 
images we use here; the others are either outside the field 
or resolved away by the much higher MERLIN resolution. 
We used SExtractor in ASSOC mode to produce an output 
catalogue of these sources. 

In order to produce a larger, fainter catalogue we 
use SExtractor in the standard mode to detect islands 
of flux above a given threshold. This requires a reliable 
noise map, which is made by estimating the local back- 
ground noise at each me sh point of a grid across the image 
(jBertin fe Arnoutdll996h . The mesh size is an important in- 
put parameter; if chosen to be too small the background 
estimation is affected by the presence of real objects; if cho- 
sen to be too large the small scale variations in the back- 
ground cannot be reproduced. We adopt a size of 32 pixels, 
which corresponds to a 2"scale. We find 3.45^iJy beam - 
rms noise in t he re sulting map, in close agreement with 
iMuxlow et al. (2005). Using this noise map, sources with 
a total flux greater than 10/^Jy (~3cr) were extracted. This 
catalogue contains 691 objects and we refer to as the 'silver' 
set of objects. 



3.2 Optical Data 

The optical data we use forms part of the Great Observa- 
tories Origins Deep Survey (GOODS) based on multiband 
HST imaging of the HDF and CDF. The HDF-N has been 
imaged in the ACS F435W, F606W, F814W and F850LP 
bands (B, V, i and z respectively) although for the purposes 
of this work we make use of only the z-band image and cat- 
alogue. 

The z— band images were observed in 5 epochs sep- 
arated by 40-50 days. In the odd numbered epochs each 
10' x 16' field was tiled with a grid of 3 x 5 individual 
ACS pointings. In the even numbered epochs the field was 
rotated by 45° and tiled with 16 separate pointings. The 
z— band image exposure time was typically 2100s divided 
into 4 exposures to ensure good cosmic ray rejection. In 
each exposure the telescope field of view was shifted by a 
small amount to allow optimal sampling of the PSF. The 
multiple epochs were then combined into a single mosaic. 
The obs ervations and image reduction are described in 
detail in iGiavalisco et all (|2004l ). We use the publicly avail- 
able SExtractor configuration files (specified for each band 
at Ihttp: / / archive.stsci.edu / pub /hlsp / goods /catalog _r2/h 
to make our catalogues; these configuration files have been 
fine tuned to minimise the number of false detections. 



4 SHEAR MEASUREMENT METHODOLOGY 



4.1 Shapelets 

Here we sum marize the s hapelet s method which is de s cribed 
more fully in iRefregierl (1200381) , iRefregier fc Baconl (|2003l ) 
and IMassev fc Refregierl ( 2005 ). In this approach, in the 
Cartesian formalism, the surface brightness /(x) of a galaxy 
is decomposed into a series of localised orthonormal basis 
functions B nit „ 2 called shapelets: 



/(*) = ^2 Ui,n 2 B ni>n2 (x-f3), 



where 



B„ 
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and where -H r m (ry) is the Hermite polynomial of order m, with 
the characteristic scale of the basis described by (3. The series 
converges most quickly if the characteristic scale j3 is chosen 
to be similar to the size of the galaxy, and the centroid of 
the object is located accurately. The sum of n\ and 712 is 
referred to as the order of the basis functions. In practice any 
decomposition has to be truncated at some order n max such 
that the decomposition yields a sufficiently accurate model 
of the galaxy while also being computationally efficient, as 
the computation time of each object's decomposition is oc 
n max- From orthonormality, we can find shapelet coefficients 
for a galaxy by calculating 



fn 1 ,n 2 — / f{ x )B ni 



(x; j3) d 2 x. 



(12) 



We use the publicly available sh a pelets software packag^3 
described in IMassev fc Refregierl (|2005l ) in order to make 
shapelet decompositions for all o ur objects. This cod e is 
well tested using op tical data (c.f. iHevmans et al.|[2006l and 
IMassev et al.ll2007l ). and we seek to extend its applicability 
to radio data here. 

The code usually fits convolved shapelet coefficients to 
a galaxy while also optimizing centroid x c , (3 and n max us- 
ing a non-linear algorithm. For the radio objects we found 
that this led to a large number of failures, due to the incor- 
rect estimation of j3 or the centroid wandering off the edge 
of the postage stamp; in order to stabilise the behaviour, 
we fix the centroid position x c to the SExtractor detection 
centroid. We also have the freedom to fix j3 to 0.4 times 
the SExtractor FWHM which we find consistently leads to 
models with reasonably low n max . Figure [1] shows some of 
the radio objects and their resulting shapelet models (still 
convolved with the beam). 



4.2 PSF Deconvolution 

To deconvolve the beam/PSF from the radio/optical data, 
we require a shapelet model of the relevant kernel. In the ra- 
dio image the restoring beam is exactly known, so the decon- 
volution process is relatively straightforward. As mentioned 
above, the restoring beam is a 0.4 arcsec circular Gaussian. 



IMassev fc Refregierl l|2005h describe the ensuing deconvolu- 
tion step in detail; briefly, the shapelets are convolved with 



In this section we describe the methods used to make esti- 
mators of the shear for all of our radio and optical sources. 



4 http:/ /www. astro. caltech.edu/~rjm/shapelets/code/ 
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Figure 1. Examples of radio sources in the HDF-N (left of each column) and their convolved shapelet reconstructions (right of each 
column). These decompositions have been carried out with fixed SExtractor centroid and f3 proportional to the FWHM. 



the PSF model, and the resulting functions are least-squares 
fit to the data. The coefficients of the fit correspond to the 
deconvolved model. 

To estimate the ACS PSF, we also use SExtractor to 
create a catalogue of stars using the SE x tracto r star/galaxy 
classifier index (see iBertin fc Arnoutsl (|l996t ) for details). 
We decompose each stellar image into shapelet coefficients; 
for each shapelet mode, the mean coe fficient is used for the 
PSF m odel. Again, the methodology of lMassev fc Refregierl 
(2005|) is used to deconvolve all galaxy images. 



4.3 Shear Catalogue 

Before arriving at a final shear catalogue, various necessary 
cuts were applied to the datasets. Galaxies with failures in 
the shapelet modelling, due to poor fits, have been re- 
moved. In the radio case we also removed all objects that 
were not resolved, i.e. FWHM< 0.4", and also applied a flux 
cut of Si. 4 > 54/i Jy to remove low level noise peaks. In the 
optical case we applied a magnitude cut of m z < 25 in or- 
der to only work with objects with S/N> 6; this cut also 
removed all unresolved objects. 

We now need to combine the shapelet coefficients of 
each object to estimate the weak shear they have experi- 
enced. We use the simp l e Gau ssian- weighted shear estimator 
given bv lMassev et all (|2007l ). 



7 = 



v/2f 2 



( fo,o — £4,0} 



(13) 



Note that the average on the denominator is taken over the 
objects once the cuts described above have been made. In 
this fashion we calculate a two-component shear estimator 



for each useable galaxy in the survey. We are now ready 
to examine the properties of weak lensing in the radio at 
current flux limits. 



5 WEAK LENSING RESULTS 
5.1 Useable number density 

In our final shear catalogue constructed as described above, 
we obtain number densities of n = 0.75 arcmin" 2 for the 
gold radio dataset, n = 3.76 arcmin -2 for the silver radio 
dataset, and n = 40.66 arcmin -2 for the optical data. It can 
immediately be seen that currently, radio number densities 
are substantially lower than those available at optical wave- 
lengths. However, it s hould also be noted how much n has in- 
creased in relation to IChang et all (|2004h . where there were 
only ~ 20 objects per square degree. With the imminent ar- 
rival of e-MERLIN and LOFAR, radio number densities will 
begin to compare well with optical number densities; this is 
crucial for weak lensing 2-point statistics, where the noise is 
inversely proportional to the number density. 



5.2 Shear Statistics 

In addition, the noise on weak lensing 2-point statistics is 
proportional to the shear estimator variance, ct 2 , so this is 
an important item to quantify, together with the overall dis- 
tribution of shear estimators. A histogram showing the dis- 
tribution of the shear estimators for the gold and silver radio 
sets are shown in Figure [2] together with the distribution of 
the optical shear estimators. We see that the bin size is nec- 
essarily large for the gold set, on account of the small num- 
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Figure 2. Normalised distribution of shear estimators for the de- 
convolved gold (solid line) and silver (dashed lines) radio objects 
and the optical (dot-dashed line) objects. 

ber of objects; nevertheless, we can immediately see that the 
variance of both radio sets and the optical set are compara- 
ble. Focussing on the silver and optical sets where there is 
more detail, we can see qualitatively that the distributions 
are similar in shape, with a peak and wings deviating from 
Gaussian form. 

More quantitatively, for the radio samples we find that 
<r 71 = 0.29 ± 0.03, cr 72 = 0.30 ± 0.03 for the gold set, and 
<t 71 = 0.35±0.02, cr 72 = 0.41±0.02 for the silver set. This can 
be compared with cr 71 = 0.326 ± 0.004, ct 72 = 0.328 ± 0.004 
for the optical objects. 

This is encouraging; it shows that radio shear estimators 
are not much more noisy than optical shear estimators at 
this flux limit. The combined effects of intrinsic ellipticity 
and measurement noise give comparable shear variance in 
both parts of the spectrum. 

In Figure [3] we show how the quantity <j 71 varies with 
flux in the radio case (top panel) and with the z-band mag- 
nitude in the optical case (lower panel). For both the gold 
and silver objects, we find that the fainter objects have a 
larger scatter compared to the brighter ones, due to the re- 
duced signal-to-noise, and the greater noise associated with 
deconvolution for small objects. The optical objects show 
a less pronounced trend, with bright objects having a little 
less scatter than the fainter ones. 

5.3 Systematic Errors 

We now wish to quantify the level of systematics present in 
the processed radio shear data. This paper seeks to assess 
the level of the problem, and to take a straightforward first 
step towards its amelioration in section [6] 

Firstly, we calculate the average shear for the entire cat- 
alogues. For the gold radio set, we find (71) = 0.035 ± 0.040 
and (72) = —0.007 ± 0.040, consistent with no overall sys- 
tematic offset. For the silver set, we find (71) = 0.072±0.021 
and (72) = 0.015 ± 0.025, indicating a systematic effect 
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Figure 3. Top panel: Shear estimator dispersion <r 71 for the de- 
convolved radio objects as a function of flux. Lower panel: cr 71 for 
optical data, as a function of z-band magnitude (with axis flipped 
for comparison with top plot). 



afflicting the estimators at this lower signal-to-noise. For 
the optical estimators we find (71) = —0.014 ± 0.006 and 
(72) = —0.003 ± 0.006, again showing an uncorrected small 
systematic due to the basic PSF correction, and/or a cosmic 
shear signal. 

As a more detailed test, we median averaged our shear 
estimates in 1 arcminute bins in x and y position on the im- 
age (Figure!?}. F° r the gold radio objects (top panel) we see 
that 71 and 72 are consistent with zero in the x direction, 
with the exception of one bin for 72. In the y direction the 
results are again consistent with no systematic. For the silver 
set (middle panel) we find 71 to be contaminated in both the 
x and y direction, with 72 showing no evidence of systemat- 
ics. For the optical objects (bottom panel) we see that there 
is a per cent level contamination in 71, with no evidence of 
systematics in 72. We also show in Figure [5] how the tan- 
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Figure 4. Median shear estimators in x and y position bins. The 
top, middel and bottom panels show results for the gold radio, 
silver radio and optical objects respectively. 



gential component of the shear 7t = 71 cos(2f?) — 72 sin(2r?) 
varies with radial distance from the centre of the image. 
We see that in both the radio cases as well as the optical 
case this is consistent with zero even on the outskirts of the 
images. 

In order to isolate the source of the systematics present 
in the silver and optical data sets, we examine the median 
average shear in radio flux bins or optical z-magnitude bins 
(Figure [S]) and also as a function of object size (Figure 0. 

Firstly considering Figure [6] we see that for the gold ra- 
dio objects (7} is consistent with zero as a function of flux. 
In the silver radio case we see that the fainter objects exhibit 
non-zero systematic shear. This may be due to anisotropic 
noise in the radio image which has not been removed by our 
deconvolution process. In the optical case we see some vari- 
ation of (7) with magnitude, with the strongest deviation 
from zero seen at the fainter magnitudes. This can be un- 
derstood as uncorrected PSF systematics due to our basic 
correction method. 

When we plot (7) as a function of FWHM (Figure[7]), we 
again see that there is deviation from zero shear in the radio 
objects. At large FWHM, a few objects with large ellipticity 
dominate the signal. More troublesome is the non-zero shear 
for the silver at small FWHM; this is a further manifestation 
of the possible anisotropic noise problem raised above. The 
optical case is mainly consistent with zero. 

A final test makes use of the correlation function £tr(0) 
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Figure 5. Tangential shear as a function of radial distance from 
the centre of the image. 



that should be zero in the absence of systematics. In Figure 
[8] we show our gold and silver radio £t r (#) correlation func- 
tions, along with the optical data and optical-radio cross- 
correlation. 

We see that, while in the optical and gold cases, the 
results are consistent with zero, in the silver case there is 
clear evidence of systematic contamination on scales ^ 2' 
and 7'. This can be fully removed by using the optical- 
radio cross-correlation, as will be discussed further in H6.2I 
in this cross-correlation case, the results are again consistent 
with zero as shown on the plot. 

In summary, then, we see that there are substantial 
residual systematics in our silver set which are not evident 
in the gold set; these may be due to anisotropic noise, and 
pose a challenge which radio weak lensing studies need to 
solve. However, systematics of this size do not defeat the 
analysis of this paper, as we shall now see. 

5.4 Cosmic Shear Constraints 

We are now in a position to use our shear estimators to 
constrain cosmological information, using shear correlation 
functions as described in ij2] recall that the correlation func- 
tions are simply related to the shear power spectrum and 
hence the matter power spectrum; equally, they will respond 
to any systematics present. 

We calculate the £+ and correlation functions for 
our gold, silver and optical sets, and their cross-correlations, 
using equation (|6]), then convert to i?-mode and B-mode 
correlation functions using equations (|SJ) and © . Our results 
are shown in Figures [9] and II 01 

In each case we see B-mode contamination on scales 
< 2"; however, on larger scales there is no evidence of 
systematics. The gold and silver radio correlation functions 
have large error bars due to the low number density; in 
these cases we have not attempted to fit cosmological pre- 
dictions. The optical-radio cross-correlation functions have 
significantly smaller error bars and it is these that we use to 
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Figure 6. (7) in radio flux bins for the gold and silver radio 
objects (top and middle panels) and in 2-band magnitude bins 
for the optical objects (lower panel). 



Figure 7. (7) in FWHM bins for the gold, silver and optical 
objects, for 71 (top) and 72 (bottom). 

constrain cosmology; however, even here we are only obtain- 
ing an upper bound on the cosmic shear signal, and hence 
on cosmological parameters. 

We calculate correlation function predictions for a flat 
ACDM model with Ho = 72 km s~ 1 Mpc~ 1 with two varying 
parameters, fi m and erg, using equations (0 to ©. These are 
X 2 fit to the measured correlation functions, with resulting 
constraints shown in Figure [TT] We find upper limits on the 
cosmological parameters, 

°*(^)°' 5 ^ 6 < 1 - 1 ^ ( 14 ) 

for the radio-gold optical cross-correlation, while for the 
radio-silver optical cross-correlation we find 

^b(^)°' 5 ^ 6 <0.8 (la) (15) 
These results do not include the cosmic variance due to the 
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Figure 8. Top panel: §t r correlation functions for the gold ra- 
dio objects overlaid with the optical and the gold-optical cross 
correlation functions. Lower panel: same plot for the silver set, 
together with optical and optical-silver cross-correlation. 



fact that the HDF-N is a small field, but show that radio- 
optical weak lensing cross-correlations are already capable 
of constraining cosmology at an interesting level. 

In these constraints we have included the median red- 
shift as one of the parameters, as we do not have a complete 
redshift sample for our radio catalogue. However, available 
data suggest that the median re dshift is near z m ~ 1 : using 
2600 spectroscopic redshifts from lBarger et all (|2008h in the 
GOODS North with the same median magnitude as our op- 
tical sample, we find a median redshift for the optical sample 
to be z m = 0.84. This is in close agreement with the median 
redshift z m = 0.85 found from NEE0 for 84 of our matched 
radio-optical objects; this can be compared with th e pre- 
dicted median redshift from the lWilman et all (|2008l ) radio 

5 http:/ /nedwww. ipac.caltech.edu/ 
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Figure 10. Cross correlated E- and B- modes, with gold-optical 
on the top panel and silver-optical on the bottom panel. 



simulations at the relevant flux threshold, which is found to 
be z m = 1.1 for Si. 4 > 50/^Jy. 



Figure 11. Cosmological constraints derived from the gold- 
optical (top panel) and silver-optical (bottom panel) cross cor- 
relation functions at z m = 1.0. The solid and dashed lines are the 
lcr and 2<r contours respectively. 

while a weak correlation provides a degree of independence 
of radio and optical lensing measurements, as we will discuss. 



6 OPTICAL-RADIO CROSS-CORRELATIONS 

We have seen above that, given the residual systematics in 
our radio shear data, it is useful to cross-correlate shape 
information between radio and optical data. The question 
then naturally arises: are radio and optical shear estimators 
well correlated for individual galaxies? 

The shear estimator of a single object is dominated by 
its intrinsic ellipticity, so radio and optical shear estimators 
will be strongly correlated if the radio and optical emission 
is aligned and has a similar ellipticity. Either a strong or 
weak correlation would be of interest; a strong correlation 
would afford a useful check for shear measurement methods, 



6.1 Observed Correlation Strength 

We matched our radio and optical catalogues, finding 123 
objects with emission in both wavelengths. In our sample 
there were 4 instances in which multiple radio sources were 
associated with a single optical one. These have been re- 
moved to avoid any ambiguity about which to compare with 
the optical estimator. We show the shear estimators for these 
galaxies in both parts of the spectrum in Figure [12] We find 
that they are in fact not strongly correlated, with a Pear- 
son's correlation coefficient of 0.097 ± 0.090. 

This result is initially rather surprising, but is confirmed 
by visual inspection of the radio and optical shapes of the 
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Figure 12. Comparison of radio and shear estimators for 
matched objects. 



objects: as an example, in Figure [T3l we show the brightest 
objects in our radio sample that have optical counterparts. 
We have overlaid the 2, 3 and 4<r contours (where a = 3.3/iJy 
beam -1 is the rms noise of the radio image) on the z-band 
optical image. Note that there is little evidence of tracing or 
alignment between radio and optical emission. This is also 
confirmed by a statistical analysis in SJ7] We will now discuss 
the implications of this result, both for weak lensing, and in 
comparison with other work in the field. 



6.2 Weak Lensing Implications 

The lack of correlation in radio and optical shapes is of real 
use for weak lensing studies. The basis of many studies in the 
field is the measurement of 2-point statistics after measuring 
shear estimators; in the weak lensing regime, we can write 
the shear estimator as 



7 = 7 + 7 +7 



(16) 



where 7 is the true gravitational shear signal, 7 1 is propor- 
tional to the intrinsic ellipticity of the galaxy, and 7 s repre- 
sents any shear contributions induced by systematic effects. 
If we then write out a correlation function, we obtain the 
terms 



(77} 



(77) + ill) + 2(7*7) + <7 S 7 S )- 



(17) 



where we have neglected terms (77 s ) and (7' 7 s ) as evidently 
zero. 

The first term on the right hand side of equation (JT7J) is 
the true cosmic shear signal that we desire to measure. The 
second term arises when physicall y close galaxies are i ntrin- 
sically aligned; this is the 'IF term lMackev et al.1 (|2002T l. The 
third term can arise when a foreground gravitational poten- 
tial tidally distorts a neighbouring galaxy, while also lensing 
a back ground galaxy; this is the 'GF term lHirata fc Seliakl 
|2004 ). The final term involves only systematic shear. In gen- 
eral, this can be large and requires a good understanding of 
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Figure 13. All matched objects with optical m z < 21.5 and radio 
flux > 140^tJy. The z-band optical images are overlaid with the 
2, 3 and 4<r radio contours. 



the experiment and sophisticated techniques to remove any 
spurious signal. 

The disentangling of the pure (77) term from the above 
equation is therefore very challenging. However, let us now 
suppose we have a joint radio and optical survey, containing 
a set of galaxies each of which has a radio or optical shear 
estimator (•> and 7, respectively), or frequently both. In 
this case, a cross correlation function can be measured: 



= <77) + (7o7> + (7r7> + (7o7r) + (7o7r)- 



(18) 



This leads to several advantages. Firstly, it should be noted 
that if the covariance between shapes in optical and radio 
is negligible, the error upon the first term is substantially 
reduced due to an effective increase in galaxy number den- 
sity, as shown by l)jarvis fc Jai 11 2008) ; we can effectively 
consider the optical and radio sources as independent galax- 
ies. This is an advantage over cross-correlating two optical 
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bands, where a strong shape correlation is fo und and little 
incre ase in effective number density results (|jarvis fc Jainl 
l2008h . 

We see that besides the cosmic shear signal, we now 
have two GI terms, one from the optical and one from 
the radio data; these will still need to be measured , 
marginalised over or nulled (|joachimi fc Schneider! l2008h . 
However, the other terms are reduced in amplitude by the 
cross-correlation: in i )6.f I we have shown that (707,*} is small; 
there is little correlation between optical and radio shapes. 

In addition, the last term will be small; the system- 
atics associated with a given optical and radio telescope 
and ensuing data reduction are so distinct that they could 
hardly be correlated . This is related to the approach of 
( ,J arvis fc Jainl [200i ) who advocate cross-correlating shear 
measurements in different exposures; however, this might 
still leave chronic systematics uncorrected, which would be 
removed by a radio-optical correlation. 

We conclude then that cross-correlating future large ra- 
dio and optical datasets could be a powerful method for lens- 
ing studies, substantially reducing the issues associated with 
intrinsic alignments and systematics. Our results in ij5]4] al- 
ready give some evidence of this but better radio data will 
be required for the cross-correlation to be fully studied and 
utilised. 



6.3 Alignment of Optical and Radio Emission 

Recently iBattve fc Browne! (|2009l ') studied the orientation 
of radio and optical galaxies in the FIRST and SDSS sur- 
veys. They found a clear excess of galaxies in which the 
major axes of the radio and optical emission were aligned. 
This may appear to be in tension with what we have shown 
in Figure 1121 but the comparison is not direct; since this 
figure shows shear estimates, it mixes orientation and flat- 
tening information. In o rder to compare more directly with 
IBattve fc Browne! (|2009h . we calculate for all matched ob- 
jects (solid histogram) the angle between radio and optical 
major axes; our results are shown in Figure fH] We also plot 
in Figure [14] the corresponding histogram for the third of 
the matched objects with the highest signal-to-noise (dashed 
histogram) . 

We see that there is only a lcr excess of the total objects 
with a small angle between radio and optical major axes. 
We also see that the histogram for the high signal-to-noise 
objects does not show any evidence of such an alignment. 
Thus there is some level of disparity between our work and 
IBattve fc Browne! (|2009l ). We not e that our sample i s muc h 
smaller in number than that of IBattve fc Browne! (|2009h . 
and as such, we are not able to usefully bin our sample more 
finely than shown in Figure IT41 

Equally it is important to note that this study involves 
highly resolved imaging using sub-arcsecond angular reso- 
lution radio and optical data. In comparison, the study by 
IBattve fc Browne! (|2009l ) uses much lower resolution data; 
they use VLA FIRST radio data with an angular resolution 
of 5 arcsec correlated with SDSS which has typical seeing 
in the range 1 to 1.5 arcsec. At these resolutions both the 
radio and optical data typically trace the extended smooth 
emission from the galaxy and hence its shape on scales of a 
few tens of kpc at the typical redshifts of these s ources . As 
such the correlations seen by IBattve fc Browne! (E)09T) arc 
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Figure 14. Histograms of the matched radio and optical galaxies 
in our sample binned according to the quantity Act, defined as 
the acute angle between the optical and radio major axes. The 
solid histogram is the full matched radio-optical sample while the 
dashed histogram corresponds to the third of the matched objects 
with the highest signal-to-noisc. 



not surprising. The objects studied here in the HDF-N are 
much fainter, and at the 0'.'4 angular resolution of this study 
we are typically resolving structure with linear sizes of just 
a few kpc. 

Additionally the types of galaxies involved in each of 
these two studies differ considerably, with the sources in 
the HDF-N sample being predominantly moderate to high 
redshift systems with intense youn g star-formation and large 
levels of optical extinction, and the lBattve fc Browne! (|2009f ) 
sources being lower redshift, more quiescent galaxies. 

Thus there are clear reasons why there may be differ- 
ences in these two studies. However the precise cause of this 
disparity is worthy of further investigation, especially as it 
has strong implications for the design of future combined 
optical and radio weak lensing studies. 



7 BLIND SHEAR MEASUREMENTS 

In this section we describe a supplementary a pproach to ra- 
dio sh ear measurements. We are motivated bv lMuxlow et all 
(|2007l ) who find 92 radio sources at a detection threshold of 
40/iJy within a 10' x 10' region centered on the HDF-N. 
This should be compared with the size of the HST ACS z- 
band catalogue, which contains ~ 13000 galaxies brighter 
than m z = 28.3 in the same region. Although a vast ma- 
jority of these sources are not detected individually at radio 
wavelengths, it is possible to statistically detect these very 
weak radio sources. Figure [15] shows the mean radio flux 
measured within a 0.75 arcsecond radius of the positions 
of ~ 8000 of these sources, as a function of their optical 
magnitude. Radio emission at the level of a few /xJy is sta- 
tistically detected from optical systems as faint as m z = 25; 
this is in good agree ment with the analysis conducted by 
iMuxlow et all (|2007n . The lower data points in Figure [T5\ 
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Figure 15. The mean 1.4GHz flux density detected within 0.75 
arcseconds of the positions of optical galaxies in the field, binned 
in z-band magnitude. Radio emission is statistically detected from 
optical galaxies brighter than m z = 25. Lower data points: results 
for apertures placed at random on the image. 




show mean flux when aperture positions are chosen at ran- 
dom. 

Since the surface brightness of these objects can be sta- 
tistically detected, we can also attempt to statistically quan- 
tify their ellipticities. In order to do so, we choose optical 
galaxies in a particular magnitude bin, rotate them so that 
their major axes are parallel, and create a composite op- 
tical galaxy by finding the median value of all galaxy sur- 
face brightnesses on a pixel by pixel basis. We then use the 
positions of the optical galaxies to extract postage-stamps 
from the radio image of diameter 3.75". We rotate these by 
the same amount as the relevant optical image, and median 
average these also. The resulting stacked optical and radio 
images are shown in Figure [16] 

As expected the stacked rotated optical objects are el- 
liptical in shape. If radio object ellipticities were well-aligned 
with optical ellipticities, we would see the stacked rotated ra- 
dio objects having an ellipticity on our plots; instead, they 
are rather circular in shape, suggesting a low correlation 
in alignment, in confirmation of our result in |6] Quantita- 
tively, we measure the ellipticity of the images in Figure [TBI 
We calculate the quadrupole moments of the images, 



/ d 2 x I(x)W{x) 
J d 2 xl(x) 



(19) 



where I(x) is the intensity profile of the object and W(x) is 
a weighting function, which here is a top-hat function with 
radius equal to a third of the radius of the postage stamp. 
The 2-component ellipticity can then be calculated as 



ei,2 = 



{Qn -Q 2 2,2Q 12 } 



hi + i 



(20) 



The calculated ellipticities are shown as a function of the 
z— band magnitude in Figure [TT] with errors calculated via 
jack-knifing of the object ellipticities in each magnitude bin. 
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Figure 16. Left: median stacked, rotated optical images. Right: 
corresponding median stacked, rotated radio postage stamps. 
From top to bottom: m z ^ 20, 20 < m z ^ 22, 22 < m z ^ 24, 
24 < m z < 26. 
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Figure 17. Ellipticities of the rotated optical (solid line) and 
radio (dashed line) composite images in Figure H"6l 



We see that the ellipticities of the rotated stacked radio im- 
ages are much lower than the optical stacked images; again, 
this acts as a confirmation of our results in |6] as it means 
that radio emission is not coherently oriented with respect 
to the optical emission. We see that at optical magnitudes 
greater than m z = 23 the error on the radio ellipticity is 
substantial, as the radio emission is not well characterized 
beyond this magnitude limit (see bottom panel of Figure 
I16[) . We have therefore applied a magnitude cut of m z ^ 23 
in what follows. 

The natural extension to this analysis is to form blind 
radio shear correlation functions. By using only the posi- 
tions of the optical objects we can extract postage stamps 
from the radio image and again calculate the ellipticites 
using the quadrupole approach above; we can then calcu- 
late the correlation functions for the simple shear estima- 
tor 7 ~ e/(2 — (e 2 )). The resulting correlation functions 
are shown as solid lines in Figure [THJ The errors from this 
technique are comparable with those found for optical-radio 
correlations in section [5741 resulting cosmological constraints 
are shown in Figure 1X91 amounting to 

-s(£f|) 'V<1.0 (la) (21) 

We can also make use of the optical shear estimators we have 
already calculated. We cross-correlate these optical shears 
with the blind radio ellipticities, showing the resulting cor- 
relation functions as the dashed lines in Figure [18] and cos- 
mological constraints in Figure [191 giving 

-8(^)°' 5 ^ 6 <1.4 (la) (22) 

The fact that these constraints are comparable with those 
from the more conventional radio-optical correlation func- 
tion in 15.41 will lead us to pursue both techniques in future 
work. 
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Figure 18. Blind radio shear correlation functions, showing the 
i?-mode (top) and B-modc (bottom). Solid line: using optical po- 
sitions only, with radio ellipticities; dashed line: cross-correlating 
optical shear estimators with blind radio ellipticities. 

8 CONCLUSIONS 

In this paper we have endeavoured to measure some key 
properties of weak gravitational lensing at radio wave- 
lengths. The motiva tion for this is the fact that to date, al- 
most all (apart from lChang et ai1l2004l ) weak lensing studies 
have been carried out at optical wavelengths. Forthcoming 
radio telescopes such as e-MERLIN, LOFAR and ultimately 
the SKA with their massively increased sensitivity, will be 
able to detect a vast number of radio sources that will make 
weak lensing accessible at radio wavelengths, so this work 
makes some early steps towards the full exploitation of those 
future surveys. 

In this work we have taken some of the most sensitive 
radio data currently available and measured the shear es- 
timator distribution of the detected radio sources. Making 
use of the publicly available optical HST ACS data we were 



Radio Weak Lensing with VLA and MERLIN 15 



2.0 




0.2 0.4 0.6 0.8 1.0 




0.5 - 



0.2 0.4 0.6 0.8 1.0 

n m 

Figure 19. Cosmological constraints from the blind-blind (top 
panel) and blind-optical cross correlation functions with z m = 
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able to directly compare the properties of the radio sources 
and the correspond ing optical galax ies. We have used the 
shapelets method of lRefregierJ (|2003al ) to measure the shapes 
of our sources while deconvolving the beam/PSF. For the 
radio data the beam was exactly known, whereas for the 
optical data the PSF was estimated using the stars in the 
image. We found a useable radio number density for shear 
analysis of n = 3.8 arcmin -2 compared to that of the optical 
n — 40.7 arcmin -2 . The radio number density is therefore 
still rather low c ompared to optical, but is much larger than 
that available to lChang et alj (|2004h ■ 

We used a simple shear estimator for our radio and op- 
tical objects and found comparable shear distributions; the 
dispersion cr 7 governs the error on the derived cosmic shear 
correlation functions, and we find its value in the optical 
and radio is comparable (<r 7 = 0.3 to 0.4 depending on the 
particular sample). 



We present various tests of systematics, including mea- 
suring mean shear as a function of image position, flux, and 
FWHM, finding little evidence of systematics in our gold ra- 
dio set, but substantial remaining systematics in the silver 
radio and optical sets, particularly at a low flux limit. In the 
silver radio case, this may be due to the effect of anisotropic 
noise affecting the measured ellipticity of objects. 

We measure shear correlation functions for the gold ra- 
dio, silver radio and optical datasets, finding B-mode con- 
tamination on scales < 2'; the uncertainty on the radio 
correlation functions is too large to provide significant cos- 
mological constraints, whereas the cross-correlation of ra- 
dio and optical shear estimators provides an upper bound 
(T8(fim/0.25) ' 5 2J; 6 < 0.8 (la, not including cosmic vari- 
ance). This cross-correlation also removes the observed sys- 
tematic error in the £tr correlation function. 

When we compared the shear estimators in optical and 
radio for the matched objects, we found very little correla- 
tion between the shear estimators of the radio source and 
the optical counterparts, with a Pearson's correlation coef- 
fic ient of ap = 0.097 ± 0.0 90. This is in apparent contrast 
to lBattve fc Browne! l|2009h who recently showed that there 
is an excess of galaxies in which the major axes of the opti- 
cal and radio emission are aligned when comparing optical 
galaxies from SDSS and radio sources from FIRST. How- 
ever, we have argued that our very different regime of depth 
and resolution accounts for this disparity. 

We have discussed how the low correlation can be ben- 
eficial in extracting a cleaner cosmic shear signal given the 
availability of both radio and optical data sets, by measur- 
ing the radio-optical shear cross-correlation function; this is 
borne out by the improvement in systematics this provides 
in our measured cosmic shear constraints. Additionally, the 
fact that the galaxies have uncorrelated shapes in the radio 
and optical data increases the effective so u rce de nsity. 

Finally, motivated by iMuxlow et al.l (|2007l ) , we made 
blind measurements of radio shear at the k nown positions 
of opt ical galaxies. We confirmed the result of lMuxlow et al.l 
|2007l ) that radio emission down to a few /^Jy is statistically 
detected from optical galaxies as faint as m z — 25. Build- 
ing on this, we proceeded to measure the average shape of 
these statistically detected objects. By stacking optical im- 
ages that had been rotated by their orientation angle along- 
side postage stamps of the radio image with the same optical 
position and rotation, we found very little association be- 
tween the ellipticities of the resulting optical stacks and the 
radio stacks, confirming our earlier optical-shear correlation 
result. 

We constrained the shear correlation functions in this 
blind approach; using either the radio ellipticities found at 
the optical positions, or using the cross-correlation between 
optical shear estimators and radio blind ellipticities, we ob- 
tain errors comparable with the more conventional approach 
above. 

The principal goal of this paper was to determine 
whether weak lensing at radio wavelengths is possible and 
promising. The improved number densities and suitable a 7 
reported here, together with the reasonable cosmological 
constraints measured, argue towards this being the case. 
Indeed, the precise knowledge of the beam for the decon- 
volution step provides radio weak lensing with a potential 
advantage over optical studies. However, a careful analysis 
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of systematics will be required, as we have shown that they 
are currently dominant unless one cross-correlates with op- 
tical data. Perhaps this is a satisfying conclusion: that in the 
future, we may not favour optical or radio shear surveys, but 
both together. 
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